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1. INTRODUCTION 



The purpose of this paper is to investigate and recommend suitable 
types of measurements on the lunar surface which will result in the 
establishment of a selenodetic control network. In this paper, we will 
define a selenocentric celestial cooi’dinate system, in which we may 
make measurements, and in turn determine an astronomic position on 
the lunar surface. Combining this astronomic position with a knowledge 
of the deflection of the vertical at this position, we could then determine 
the coordinate of this position in a selenographic system. 

Therefore this paper will be divided into two primary areas which 

are: 

(1) Establishment of a selenocentric celestial coordinate system 
and measurements in that system from the lunar surface. 

(2) Determination of the deflections of the vertical on the lunar 
surface. 



v 
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2. ESTABLISHMENT OF AND MEASUREMENTS IN THE LUNAR 
ASTRONOMICAL COORDINATE SYSTEMS 

On the earth most of the geodetic astronomic calculations are made 
using a geoequatorial system where the star coordinates are right ascen- 
sion and declination. There are a great number of catalogs, maps and 
almanacs giving the coordinates of the stars and other celestial bodies 
in this system. If the coordinates in this system could be easily trans- 
formed to a lunar equatorial system, all the methods of geodetic astronomy 
could be employed on the moon for selenodetic purposes. In this section 
we will show such a transformation complete with error analysis, while 
also discussing such related factors as time, location of lunar pole stars 
and observations in the horizon system. 

We will also show some methods which may be used on the lunar 
surface to determine astronomic position. These methods will be dis- 
cussed in relation to the instrmnent requirements, and accuracy of the 
system as a whole. Included among these methods is an inertial system, 
which is presently installed on the lunar module, and which in itself can 
establish a position on the lunar surface. 

The methods and instruments presented are not assumed to be a 
complete listing, but are based on practical factors which may obviously 
limit other methods because of physical size or complexity of equipment, 
time required to achieve an astronomic position, or excessive burden 
being placed on the astronauts, either during training or during the 
operation itself. In this light we would most like to find a method which 
is as conservative of the astronauts time as possible (in both training and 
employment), achieves satisfactory results, and does not involve undue 
calculations or recording procedures. In conjunction, we wish the 
instruments involved in this method to be simple, portable, and of 
course accurate. 
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2. 1 Lunar Astronomy 

Kolaczek [14] defines a selenoequatorial cooi'dinate system (lunar 
right ascension system) where the axis of rotation of the moon defines 
the celestial pole, which is the basic direction, and the celestial 
equator, which is a plane perpendicular to the rotation axis, containing 
the center. The latter is defined as the basic plane. This system is 
affected by precession and nutation (physical libration) of the lunar axis 
of rotation. The selenoequatorial system is analogous to the geoequatorial 
system on earth, but the selenoequatorial coordinates of the stars change 
their values faster than the geoequatorial coordinates due to the lunar 
precession being 1,360 times faster than the earth’s. The selenoequatorial 

coordinates are: 

O 

a. - Lunar right ascension measured on the lunar celestial 
equator from its ascending node on the ecliptic. It is 
measured to the east from 0° to 360°. 

d - Lunar declination which is the angular distance from the 
lunar celestial equator. It is measured from 0° to 90°, 
positive in the northern celestial hemisphere, and negative 
in the southern. 

Now we will take a coordinate system with known coordinates and 
transform these coordinates into the selenoequatorial coordinate system. 

The translation of any coordinate system from one point in space to 
another changes the values of the spherical coordinates of fixed points 
on the celestial sphere. This change is caused by the parallax or the 
translation of the origin of a coordinate system and the aberration caused 
by the different motion of this newly translated system. The translation 
of mean earth oriented coordinates into mean lunar oriented coordinates 
does not require consideration of parallax if we consider only the coordi- 
nates of stars, as we may omit the parallax of the earth-moon distance, 
due to the great distance of the stars, as compared to the earth-moon 
distance. The transformation of the apparent coordinates however requires 
consideration of the influence of the moon’s aberrations. 
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2.1.1 Transformation of Mean Geoequatorial into Mean Selenoequatorial 
Coordinates 

The transformation of mean geoequatorial into mean selenoequatorial 
coordinates can be performed in two steps [14]. First the mean geoequa- 
torial coordinates O' , 6 are transformed into the mean ecliptic coordinates 
X, /3 by the well known formulas: 

sin/3 = cos € sin5 - sin € cos 6 sino: 
cos/3 cosX = cos <5 cos a 
cos$ sinX = sin€ sin 6 + cos € cos 6 sin# 
where € the obliquity of the ecliptic is 

€ = 23°. 452294 - 0°. 0130125 T - 0°. 00000164 T 2 + 0°. 000000503 T* 
where T denotes the time measured in Julian centuries of 36,525 ephemeris 
days from the epoch 1900 January 0.5ET = J. D. 2415020.0. 

Next the transformation formulas of mean ecliptic coordinates into mean 
selenoequatorial coordinates. 

sin dm = sin/3 cos I + cos/3 sinl sin ( X - £2) | 

cos dm sin oc ? = sin /3 sin I - cos /3 cos I sin ( X - £2 ) r 2.1 

cos dm cos = - cos/3 cos ( X - £2) 

where I is inclination of the lunar equator to the ecliptic 

I = 1° 32’ 01" 

/ 

£2 The mean longitude of the ascending node of the lunar orbit on 
the ecliptic is 

£2- 259°. 18305556 - 1934°. 14200833 T + 0°. 002078 T 3 + 0°. 000002 T 3 
where T is defined as above. 

2.1.2 Transformation of Mean Selenoequatorial Into Apparent Selenoequatorial 
Coordinates 

The transformation of mean into apparent coordinates is given in 
a form [14] similar to that used for the earth’s equatorial coordinates 
£2 £2 

^ sapp = o; m + Ai ai + Bi bi + Ci ci + Di d^ + T a u a 
d sapp = dm + Ai al + Bi b] + Ci cl + Di dl + T a y d 
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Similarly to Bessel’s equations on the earth, Ai and Bi are the effects 



of precession and nutation to date from the determination of the mean 
coordinates. The Ci and Di terms are the effect of lunar annual abberation. 
Ther a u ( d/ a ) terms are the corrections for proper motion. The 
ai,bi,ci,di,ai ,bi ,cl ,dl coefficients are the lunar Besselian star constants. 
Specifically: 

Ai = - (Pst + CT ) sin I 
Bi = p 

Ci — - K sin (Lap - Ci+ 180°) cos I 
Di = - K cos (Lap - 0+ 180°) 



where: 



Po 



t 



I 






K, k 



L s ,li 
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La 



P 



- Daily motion of the ascending node 
of the moon's orbit on the ecliptic 
( = - 0°. 0529539222) 

- Number of ephermeris days from the 
epoch to of the mean coordinates o: Q , d a 

- Inclination of the moon's equator to the 
ecliptic ( =1° 32' 01") 

- Physical librations in inclination and in 
ascending node respectively (see Appendix #3) 
Constant of the earth's and of the lunar 
aberration respectively 

K = k ( 1 - Vi cos ( L s - li)/V e ) 

(k = 20". 4958) 

Longitude of the sun, the moon and the 
moon's apex in the motion around the 
sun respectively 

Lap = L s - 90° + A A 
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Ve.Vl 



Aa = tan * 



*( 



Vx sin ( L s - li ) 

V e ~ Vi cos ( L» - lx ) 



) 



- Velocity of the eal'th's and the moon's 
orbital motion, respectively 



(V e = 29.8 km/sec 
Vx = 1.02 km/sec) 



The lunar Besselian star constants are: 

a 

ax = ( cot I - sin oc B tan d B ) 

, Q, 

ax = - cos oc „ 

a 

bx = cos oc a tan d B 

r Cl 

bx = - sin oc a 

n 

ci = - cos oc 0 sec d B 
t Cl 

Cl = sin oc 0 sin d B + tan I cos d n 

, O 

di = sin a o sec d B 

t Cl 

di = cos oc a sin d B 

The proper motion will take on similar values to those computed 
on the earth: 



r a = t 

jUa, /i d = The components of proper motion of the star 
2.1.3 Lunar Horizon System 

Measurements on the lunar surface will most likely be made in the 
horizon system. In the moon's horizon system, altitude aih and azimuth Ait 
are defined in the same way as that of the earth's horizon system. The 
transformation between the lunar horizon system and the lunar selenoequatorial 
system similar to that given by Mueller [16] is as follows: 

Cl _ 

cos d sin ( LHA^-.X - oc ) = - sin Zi sin Aih ^ 

sin d = cos Zx sin <5 + sin Zx cos A cos $ >2.2 

. Cl . \ 

cos d cos ( LHA^ -X - oc ) ■- cos Zx cos $ - sin Zx cos Axh sin'SJ 
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where LHA 0 - The lunar hour angle of the node measured from the 
reference median 

<i> - The lunar astronomical latitude 
X - The lunar astronomical longitude 



Zi — 90 - aih 



The above equations can be derived with the use of Figure 1 below: 



NCP 




h = LHA^ - X - a 



Figure 1. - The Astronomical Triangle 



2. 2 Lunar Time 

As we have defined the lunar right ascension system in terms of 
ascending node, it would only be appropriate, that we define the lunar time 
system in terms of the draconic month. The draconic month is defined by 
two successive passages of the moon through the same node of its orbit. 
The draconic month is equal to 27.212220 mean solar days. 

The expression of selenographic longitude in terms of time would be 
as follows: 



1 month 


= 27.212220 


mean solar days 


- 360 


1 day 


= .9070740 


mean solar days 


= 12 


1 hour 


= .03779475 


mean solar days 


= 30 


1 minute 


= .0006299125 


mean solar days 


= 30 


1 second 


= .0000104985 


mean solar days 


= 0 
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Figure 2 below shows the relationship in draconic tunes. 




LMDT Lunar mean draconic time 

LADT Lunar apparent draconic time 

RMDT Reference mean draconic time 

RADT Reference apparent draconic time 

eq.N Equation of the node 



Figure 2. - Draconic Time System 



From Figure 2 it is evident that: 

LMDT + eq. N = LADT 
and LMDT - RMDT = X 

where X - The reduced astronomic longitude with the positive 
sign to the east 

and eq. N - The lunar equation of the node which is due to nutation 
(physical libration in longitude T and in node cr ) (see 
Appendix C). 
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From the actual rotation angles of the moon, we find that eq.N = T-a . 



We will assume in this section that the meridian is not affected by 
polar motion and that the rotation of the moon is free from short periodic 
irregularities . 

The hour angle of the node is equal to the lunar right ascension of 
an object at the instant of its upper transit over the meridian from 
which the hour angle is reckoned. Therefore, draconic time can be 
determined by observing transits of stars. 

The lunar time intervals are almost equal to the terrestrial time 
intervals of the same name (91% of the terrestrial values), which is 
still advantageous from the psychological point of view because a change 
to such a time system would be almost unnoticeable to man. The 
meridian to use for a change of dates could be the 180° meridian. 

2.3 Lunar Pole Stars 



As stated by Gurevich [5] the lunar poles describe circles in the 
celestial sphere with an angular diameter of 3° . This circle is centered 
on the ecliptic pole, with a period of 18.6 years. Therefore, first 
certain stars, then others, will be most suited for use as pole stars. 
Tables of altitudes and azimuths must be compiled for the stars closest 
to the poles at a particular time. The mean north pole 
of the moon described in the earth right ascension system would be 
approximately, a = 18 h , 6 = 66°33' . We see from the FK4 catalogue 
that the following stars from Table 1 are near the north pole. 
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Table 1. Lunar North Pole Stars 



Catalog Number 


Magnitude 


a 


6 


659 


5.21 


17 h 32“ 


Oi 

o 

o 

00 

to 


701 


6.00 


18 h 36“ 


65° 27' 


664 


4.87 


17 h 37° 


68° 46' 


685 


5.03 


18 h 13“ 


64° 23' 


and approximately 


a = 6* , 5 = 


- 66° 33' for the 


south pole in Table 2 


from the Boss catalogue. 







Table 2. Lunar South Pole Stars 



log Number 


Magnitude 


a 


6 


6927 


5.3 


5 h 32“ 


- 64° 15' 


7246 


4.5 


5 h 44“ 


- 65° 45' 


7384 


5.2 


5 h 49“ 


- 66° 55’ 


7813 


5.8 


6 h 06“ 


C\1 

0 

o 

to 

to 

1 


7946 


4.9 


6 h ll“ 


- 65° 35’ 



With a knowledge of the relationship of the position of the lunar pole to 
those stars the astronauts can make a rough north alignment with a probable 
accuracy of less than 1 minute while using only a gun sight type sighting 
device . 
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2.4 Measurements on the Lunar Surface 



Presented here are four basic types of measurements which can be 
employed on the lunar surface to determine astronomic position. 

One method is the measurement of altitudes or zenith angles of at least 
two stars separated in azimuth, combined with the known celestial coordinates 
of a star and a knowledge of time. This is a method which can be employed 
with an ordinary sextant or theodolite. 

A second method is by determination of a relative bearing and altitude of 
a star combined with information from inertial and gravitational sensors, which 
is the method used by the Alignment Optical Telescope (AOT) on the lunar 
module. 

A third method is by measurement of relative azimuth angles between at 
least three known stars. This technique can be used with an ordinary theodolite, 
but also lends itself well to the use of an automatic instrument. 

A fourth method is by star photography, which can use either the first 
or third, or other methods to determine position. Here we can use a photo 
theodolite, or other metric camera. 

In the following section these four methods will be discussed in relation 
to their ease of employment on the lunar surface, instrumental adaptability 
and accuracy, and other pertinent facts. 

2.4.1 Method 1: Measurements of Altitudes 

The measurement of altitudes or zenith angles in order to determine 
astronomical positions can be accomplished by the use of a sextant or 
theodolite. The obvious advantage of this method lies in the fact that both 
sextants and theodolites are comparatively light equipment, which could 
feasibly be included on a lunar mission. The sextant is a familiar 
instrument to the astronauts, as they already have been trained in its use. 
Although, the training required to allow the astronauts to effectively use a 
theodolite could be accomplished in a short period of tune. The big 
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advantage in using this equipment is that astronomical positions can be taken 
at several points. 

As stated in the Geonautics study [7], it is not likely though that a 
theodolite would be practicable for the initial missions. The time and effort 
required to set up and measure is far out of proportion to the many other 
tasks required of the astronauts. It would also be necessary to adapt the 
fine pointing and adjusting screws for use with the pressure gloves the 
operator must wear. In addition, precise sightings through the helmet 
face plate would be most difficult. 

A hand held sextant may be used to determine a position on the lunar 
surface. The altitude of a star is read directly from the scale. With the 
aid of sub-divisions and a vernier, altitudes can be read with a first-class 
instrument to one-tenth of a minute of arc. 

The problem which arises with the use of a simple sextant is in the 
determination of the lunar horizon. In order to eliminate this problem, 
we can consider the use of an artificial horizon or a bubble sextant, so 
that we then need only define the local vertical instead of the horizon. 

The artificial horizon is a rectangular shallow basin of mercury. The 
observer places himself so that he can see the image of the body, whose 
altitude is to be measured, reflected in the mercury. He points the 
telescope toward the artificial horizon so that he uses a second image of 
the heavenly body in place of the horizon. By moving the index arm, the 
two images are made to coincide. The circle reading is twice the altitude 
which thereby reduces the precision of the instrument to 12 seconds of arc. 

A bubble sextant (bubble octant) [17] contains its own artificial horizon. 
In measuring altitudes, the image of the celestial body is brought into 
coincidence with the image of a circular level bubble which forms the 
artificial horizon. The altitude is then read on a micrometer drum. The 
fine drum is divided into 2 minute intervals which can be interpolated to 
the nearest 1 minute. In a v bubble sextant such as the Bendix "Pioneer", 
provision is made to record a series of altitudes automatically. An 
average reading can be determined from these readings. 
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These instruments are not capable of first order geodetic accuracy. 

The positional error can be a maximum of 50 meters from the 
sextant with the artificial horizon, and 600 meters for the bubble sextant. 
This, of course, does not take into account, sextant internal errors, 
bubble errors, and errors in the star tables. In addition, the recording 
of each measurement using the common sextant will burden the astronauts. 
With both instruments each star must be identified by the astronauts, 
leading to further possibilities of error. Tune recording must also be 
done manually. 

Using this method (see Figure 1) we can find azimuth from the 

formula as given by Mueller [16] 

. . - cos d sin h 

sm A = ; — — 

sin Z 

where h = LHA^-^ - 

We can find the latitude from the basic formula as given by Mueller [16] 
cos Z = sind sin^+cosd cosh cos^ 

and we can determine longitude from the formula as given by Mueller [16] 

, cos Z - sin d sin^ 

cos h = : rr 

cos d cos 5* 



2.4.2 Method 2: Measurement with the Alignment Optical Telescope 

A relative bearing and altitude of a star can be determined by the 
Alignment Optical Telescope (AOT) which is currently installed in the Lunar 
Module (LM). An astronomic position on the lunar surface can be established 
by the use of this instrument in conjunction with the Inertial Measurement 
Unit (IMU) and the Primary Guidance Navigation System (PGNCS) accel- 
erometers. The great advantage of this method lies in the fact that no 
additional equipment or training of the astronauts is necessary to obtain an 
astronomical postion. An apparent disadvantage is the fact that only one 
astronomical position can be determined during each landing due to the 
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fact that the equipment is fixed in the lunar module which in turn is fixed 
on the lunar surface. 

The AOT as described in the Grumman manual [4] is a unity power 
periscope type device with a 60° conical field of view. Its primary pur- 
pose is to provide a method of establishing an inertial reference with 
respect to celestial bodies for inflight and lunar surface operations. It 
is operated manually by the astronauts and is mounted on the PGNCS 
navigation base to establish a mechanical alignment and common reference 
with the EMU. The AOT has a moveable shaft axis which is parallel to 
the LM "X" axis and a Line of Sight (LOS) axis which is the center of 
field of view and fixed at 45° from the LM + "X" axis. 

The AOT LOS is fixed in elevation and moveable in azimuth to six 
detent positions. These detent positions are located at 60° and are 
selected manually by turning a detent selector knob on the AOT. Three 
of the positions can be used operationally for star sightings and are called 
forward, left, and right (see Figure 3). The forward, or zero detent 
position places the LOS in the x,z plane looking forward and up. The 
right, or +60° position places the LOS 60° to the right of the x, z plane and 
the left or -60° position places the LOS 60° to the left of the x,z plane. 
Each of these positions maintain the LOS at 45° from the LM + X axis. 

The remaining three detent positions are storage positions with the 
position 180° from the forward position being the normal storage position. 
The reticle control makes it possible for manual rotation of the reticle 
for use in lunar surface alignments. The AOT reticle pattern consists of 
cross hairs and superimposed archimedes spiral (see Figure 4). 

The vertical cross hair is an orientation line designated the "Y" line 
and is parallel to the LM X axis when the reticle is at the 0° reference 
position. The horizontal cross hair is an auxiliary line designated the "X" 
line and is perpendicular to the orientation line. A one turn spiral is 
superimposed from the center of the field of view to the top of the vertical 
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60° LEFT 




TO I VIEW 



FIGURE 3. - AOT Line of Sight Positions 
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FIGURE 4 . - AOT Reticle Pattern 



cross hair. The reticle is edge illuminated. 

A counter located on the left side of the AOT, displays an angular 
readout of reticle rotation. The counter reads in degrees to within ±0.02 
degrees or ±72 seconds. Interpolation is possible to within ±0.01 degrees. 

On the lunar surface, the selected star can be sighted in either the 
AOT left, right or forward positions. The astronaut using. the manual drive 
knob, adjusts the reticle to superimpose the orientation line on the target 
star. This reticle angle displayed on the AOT counter is then manually 
inserted into the Lunar Module Guidance Computer (LGC) via the Display 
Console (DSKY). This provides the computer with the star orientation 
angle (shaft angle). The astronaut then continues to rotate the reticle until 
a point on the spiral is superimposed on the target star. This second 
angular readout is then entered into the LGC via the DSKY. The AOT 
detent position and the star code number are also inserted into the LGC 
via the DSKY. The LGC can now calculate the angular displacement of the 
star from the center of the field of view by computing the difference between 
the two counter readings. Due to the characteristics of the reticle spiral, 
this Aangle is proportional to the distance of the star from the center of 
the field of view. Using this Aangle and the proportionality equation, the 
computer is then able to calculate the trunnion angle. During the lunar 
surface phase of the lunar landing mission, the Mission Control Center 
and the LGC use the unit gravity vector as measured by the LM PGNCS 
accelerometers to align the LM platform and to estimate the latitude and 
longitude of the LM on the lunar surface. 

The AOT star alignment establishes an inertial reference by employing 
the following formulas: 



Si 




X 


(Si X S 2 )^x Si 




Y 


Si X s 2 




z 



[A] [MCI] 
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Where Si and Ss are known unit vectors of stars in moon centered 
inertial space and MCI is the unknown coordinates of the spacecraft with 
respect to the moon centered inertial system. 





u Sim _ 




X 


Sb = 


Sim X S2 M/X Si„ 




Y 




SlM X S2M 




z 



[B] [SM] 



Where Sim and S 2 M are measured unit vectors of stars in stable 
member space of the stars, and SM is the assumed stable member 
coordinates referenced to the moon centered inertial coordinate system. 

[Sa] = [Sb] 

[A] [MCI] = [B] [SM] 

[SM] = [B]" 1 [A] [MCI] 

['REFSMMAT] = [B]' 1 [A] 

Therefore: MCI = [REFSMMAT] T [SM] 

Upon comparison with the gravity vector and an estimate of the lunar 
radius at the point, a position on the lunar surface can then be estimated. 

Considering only the precision of the counter reader (with interpolation) 
we see that the positional accuracy may be' as poor as 152 meters 
or 498 feet. If we considered the errors in the determination of the unit 
gravity vector, and the errors in the star tables, we would obviously find 
that this system is not adequate for determining accurate positions on the 
lunar surface. Considering a minimum standard error of 20". 8, we will 
eliminate this system from consideration for geodetic purposes. 
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2.4.3 Method 3: 



Measurement of Relative Azimuth Angles 



Method three involves the simultaneous determination of longitude, 
latitude, and azimuth based on the motion of stars in azimuth. An 
ordinary theodolite can, of course be used, but it again suffers the 
drawbacks mentioned in Method 1. 

This method though, lends itself to a concept which can allow 
automatic determination of astronomic position and azimuth. An 
"Electronic Zenith Camera" (EZC) is currently under development by 
the Control Data Corporation of Minneapolis, Minnesota. It is described 
by Carroll as a portable sensor-computer combination which performs 
both detection of stars and calculation of astronomic position [1] [2]. 

Stars are detected by their transits across a fan of radial slits placed 
at the focal surface of a vertically oriented optical system. 

The principal advantages of this system is that once set up, the 
system can determine its astronomical position on the lunar surface 
without further attention from the astronauts. The system itself can 
be easily set up on the lunar surface. This set up includes approximate 
north orientation, and approximate levelling. If necessary, the system 
could be left to operate after the astronauts have left the lunar surface. 
Radio transmission of current updated positions could be relayed to 
earth if necessary. 

It is considered that astronomical longitude and latitude can be 
determined to within 0.3 (la) seconds and azimuth to within 3 (1CT) 
seconds (not considering star catalogue errors) which is certainly 
acceptable accuracy at this time. No instrument of this type has 
currently been manufactured, but all advance planning has been completed 
and the necessary technical skill is available to produce this instrument. 

The EZC can be broken down into the system block diagram as seen 
in Figure 5. 
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Figure 5. - EZC Block Diagram 



We will examine each component starting with the sensor 
(see Figure 6). The lens system is an f/3.3 system having a 3 

centimeter aperture and a 10 centimeter focal length. The lens system has 
a 16 degree total field of view (eight degrees maximum off-axis angle), and 
can form an image of a point source restricted to a tangential dimension 
no greater than 10 seconds of arc but not restricted by the radial dimension . 
It will be necessary to shield the lens from such bright bodies as the sun 
and earth. 

The reticle has 100 radial slits (at the focal point of the lens) about 
1.4 centimeters long and 4.8 microns wide. 
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PHOTOMULTIPLIER (2) 




Figure 6. - EZC Sensor Unit 
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Directly behind the slit reticle is placed an end-window photomultiplier 
tube which collects all light passing through the slits. The output of this 
tube is amplified, filtered, and processed so that only the maximum 
signal (peak pulse) will be recorded and timed. The time input can be 
any commercially available temperature compensated crystal oscillator. 

Looking at the control unit, (see Figure 7) we find that the heart of 
the system is a miniature general purpose digital computer from which 
the sensor is controlled, and from which the operator receives status and 
calculated information. Also included are a precision oscillator (crystal 
oscillator), a power source, and the necessary interface electronics. The 
level reference is a mercury pool auto collimator. Miniature electronic 
autocollimators can reflect a beam off the surface of a completely enclosed 
mercury pool (dust free surface), the surface of which is covered by a 
silicone fluid which damps out ripples. If the unit is to function during 
the lunar night, it will be necessary to heat the mercury pool, due to the 
lunar temperature being below the freezing point of mercury. 

The DC reversal motor is under command of the control unit. The 
motor drives the instrument about the vertical axis against one of two 
stops placed 180° apart. The reversal bearing is composed of dissimilar 
metals (brass and bronze) thereby eliminating the possibility of micro- 
scopic welds, which would inhibit accuracy. 

Prior to operation, the astronaut must level the instrument to within 
30 seconds. Within this range the automatic level reference supplies 
information to the computer so as to measure the angle between the 
instrument vertical axis- and the true local vertical. The astronaut must 
also make an initial north (south) - alignment of the instrument. This 
can be accomplished by using a gunsight mechanism on the equipment and 
sighting with the aid of near polar stars such as previously described. 

This wall yield an initial approximation of north, with the exact azimuth 
to be solved for. v 
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Figure 7. - EZC Control Unit 
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We arc interested in detecting stars from +1. to +5.3 magnitude. 

It has been determined that the best results are obtained with a star 
transiting at least one slit each minute. In order to achieve this rate 
on the lunar surface, we must either wait approximately 27 tunes longer 
for a solution than we would on earth (approximately 27 hours) or have an 
internal spinning reticle which would reduce the solution tune but would 
introduce additional errors into the system (approximately 1-arc second). 

Star identification is based on an approximate knowledge of the 
observer's position, and north alignment. The computer determines the 
expected pattern of transit times from the stored star catalog. This 
computed pattern and the measured time pattern are then superimposed, 
and these target pairs which correspond closely and uniquely in magnitude 
constitute tentative identifications. These are then used to adjust the 
assumed position and azimuth and thereby to adjust slightly the computer 
time pattern so as to obtain a better fit with the measured pattern. From 
here the final solution sequence preceeds by solving a linearized version 
of the following equation derived from equations 2.2: 

Q Q, 

cot( A 0 + /3j ) sin(LHA^ - X - a t ) -sin'&cos (LHA^ - X - a i ) 

+ tan dj cos = 0 



where 




A 0 


The azimuth of the Oth slit (dependent on north alignment) 


ft 


The fixed angle between the Jth slit and the Oth 




slit 


LHA tt ’ 


Apparent lunar reference meridian hour angle of the node 
(LADT) 


X 


Longitude 


q 


Lunar right ascension of the i th star 




Latitude 


di 


Lunar declination of the i th star 
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This linearized version is then 



0 = sin Aj cot Z d*? + cos =?( tan*? - cosAj cot Z ) dh - dA 0 



where Z is the computed zenith distance 



Cl 



h = LHA^ - X - ol j 



Aj — A 0 + fi i 



and 



cot Z 



( sin d sin<& + cos d cos h cos*?) sin A 1 



- cos d sin h 



The measured quantities are azimuth and time. Latitude and longitude 
of the station are assumed. 



Since there are more measurements than unknowns, we may use a 
least square solution, with a Kalman filtering technique [10] based on the 
Newton-Raphson iteration formula, so that matrix inversions will not be 
required. This iterative procedure is based on the following equation: 



Xn + i — Xn 



t [ aN Xn + Lm ] 

Cn a N » 

[ 1 + aw Cm a N ] 



where : 





'd*?' 




~ 0 " 


matrix Xn = 3 X 1 = 


dh 


X 0 = 


0 




_dA 0 _ 




_0 



matrix a N = ia 3 = ( sin Aj cot Zi - cos*? (tan*?- cos Aj cot Zi) - 1 ) 

Cl 

Ln = cot (A 0 + Aj ) sin (LHA^ - X - a i ) - sin*? cos(LHA^ 

- X - a. i ) + tan di cos *? 
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and matrix C N +i = C N - _ c,j aN aN Cn — 

[ 1 + aw Cm aj ] 



where matrix C 0 = 3 C 3 =(A’ A) -1 

where the rows of A are the rows of ai to as . 



We can see that the incorporation of a new measurement does not 
require a new matrix inversion. Only one inversion is required and that 
is for Co which is only an inversion of a 3 x 3 matrix. After the first 
five marks no old observation information need be saved. The values of 
the elements of matrix C are indicative of the quality of the results. 

Matrix C is the covariance matrix, but only under the assumption that 
there are no systematic errors in the residual vector L . There are 
systematic errors in the residual vector L , due to the timing error 
( r libration ) and the errors in star right ascension and declinations. 

In order to cancel alignment errors in the optical - reticle system and 
normality errors between the instrument axis, the instrument is rotated 
180° , and the astronomic position is again determined. The average of 
these positions is then used. The level status of the reversal bearing 
is also measured and it is applied as a final correction. 

2.4.4 Method 4; Measurement with Photography 

The determination of astronomical position by photography can be 
accomplished with a photo-theodolite or other metric camera using methods 
such as one or three. The drawbacks of a photo- theodolite are the same 
as for an ordinary theodolite and thus would tend to make it impractical 
for early missions. 

As described in the Geonautics study [7] , a metric camera could be 
used as a zenith camera or otherwise fixed in elevation to photograph stars. 
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Various sophisticated systems such as gravity-operated automatic leveling 
devices with gimballed camera support, electronically operated automatic 
leveling devices, gyroscopic devices, etc. , could be considered to level the 
camera. The simplest and most feasible device would seem to be the 
conventional spirit bubbles. These bubbles could be placed on a camera 
support such as a tripod, so that the camera, when placed on the support, 
could be theoretically leveled within the range of two orthogonally placed 
precise level vials inside the camera body. These precise levels can be 
photographed during each exposure, thus providing a pennanent record of 

the level reference. Drawbacks in this system tend to lie in the fact that 
the bubbles are affected by differential heating from external light, so that 
they must be protected both from radiant heat and from light. The bubble, 
of the sensitivity required in the camera, would be termed a "slow" bubble 
in a conventional sui’veying instrument. It would be even slower on the 
lunar surface due to the moon’s smaller gravitational force. 

This method provides a permanent and true record of field data as 
film can be removed from camera and returned to earth. Most advantageous 
is the fact that a great number of measurements can be taken in a minimum 
amount of time after the initial set up. There is a requirement for extra 
equipment such as a tripod, but it could be retractable and made of 
aluminum or some similar material Timing of each photo would also be 
required, but this could be done by an internal clock, which would furnish 
the time to be included in the edge data. 

2.5 Error Analysis 

Let us first consider the errors which contribute to the lunar apparent 
selenoequatorial coordinates. In the transformation of the mean geoequatorial 
coordinates to the mean ecliptic coordinates, the relatively small errors in 
Q! , 6 , and € can be considered negligible. Therefore, and A can be 
considered errorless. In the transformation of the mean ecliptic coordinates 
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to the mean selenoequatorial coordinates we consider that the error in ft 
is negligible, but there is definitely an error in inclination (I). Let us 
make a conservative estimate of 20” for this error. Differentiating 
equations (2. 1) we see that the results are similar to that given by 
Mueller [16] . 

ft 

dd a = cos a o sin I d ( X - Cl ) + ( sec dn sec P cos I - tan d D tajiP ) 

o ft 

dP + sina m dl 

but d ( X - ft ) = 0 
dj3 = 0 

ft 

therefore ddn - sina 0 dl and therefore the maximum value of 

• ft ft 

dd n = sina m dl = 20". 0 sina m 

ft 2 

and da n = (sec d n cos I - secd 3 tandn sin/3) d (X - ft ) 

ft o a ft 

- cos a E sec p sec d D sin I d p - tan d a cos a B d I 

but again d (X - ft ) = 0 

dp = 0 

ft ft ft 

therefore da B = -tand 0 cosa B dl - - 20". 0 tandn cos 0 



Considering the transformation of mean into apparent selenoequatorial 
coordinates, we will use the following equations: 



da 

+ 



+ 



ft 

sapp 




+ 



5 Ai ai 

Sct 



dcr + 



3 Ai ai 
3a ft 



da 



ft 

n + 



SBibi 

3 p 



dp + 



3Bibi 
3 a ft 




, 3CiCi 
3a ft 

HI 



ft 



da B 



+ 



, 3Didi 
3 a 




3Ai ai 

31 

3Bibi 

3 I 

3 Ci ci 

31 

d Pi d! 

31 



dl + 
dl + 
dl + 
dl + 



3 Ai ai 
3 do 

3 Bi bi 
3 d 0 

3 Ci ci 
3dn 

3 D x di 
3 d B 



ddn 

dda 

ddn 

ddn 



(2.3) 
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and a similar expression for ddsapp, considering in both cases a 
negligible error in proper motion. 

In the following differential expressions we will assume: 



cos I =1 

d D = Limited to 45° 

p = 1* 

O = 1 ' 



The differential quantities will be assigned the following conservative 
values: 



da = 


10" 


= 


.000050 


radians) 




dp 


10" 


= 


.000050 


radians/ 


(estimate from [13] ) 


dl 


20" 


= 


.000097 


radians 




ddo = 


20" 


- 


.000097 


radians 




dr 

, Cl 
dee B - 


15" 


= 


.000014 


radians 


(estimate) 


20" 


= 


.000097 


radians 





Ai ai = (Pot + cr ) (cos I - sin I since m tan d m ) 



3 Ai ax 
3a 



da 



(cos I - sin I since „ tan d D ) d a 

Cl 

( 1 02676 since B tand D )da = da 

10 ". 0 



3 Ai ai 
3 I 



dl 



(Pot + a ) (cos I since B tan d B - sinl)dl 
(-.0015059) (-. 97324)dl 
0".0 



3Ai ai 
3 do 



ddo 



(Pot + a) (sin I since B sec 2 d D )ddo 
(-.0015059) (. 05352) ddo 
0".0 
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da? 

9a u 

a 



(Pot +ct) (-sinl cos a a tand = ) da B 

Cl 

(-.0015059) (-.02676) da 0 
0".0 



Cl 

Bibi = p cosa D tandn 



3Bi bi 

9p 



dp 



= cosa B tan d 0 dp 

= 10”. 0 cos a^ tan d a 



3 Bx bi 

31 



dl 



= 0".0 



3Bi bi 
3 da 



ddm 



Cl 

= p cosa m 



3 

sec 



d n ddo 



= . 00058 dd m 

= 0”.0 



3 Bibi 
3a§ 




Q Cl 

p sin a n tan d D da m 

Cl 

- (.00029) da a 

0".0 



Ci ci 



3 Ci Ci 

3 I 



dl 



o 4 L 

= K sin ( Lap - S2 + 180 ) cos I cos a c sec d D d I 

. Q 

=-K sin (Lap - f2+ ISO ) sinl cosa - sec do dl 
= ( . 000097 )(!)(. 02676 ) ( 1 ) (1. 4) d I 



= 0”.0 



3 Ci Ci 
3 da 



dda 



_ O 4 6 

= K sin (Lap - Cl + 180 ) cos I tand : secd D cosa m dd m 
= (.000097) (1) (.02676) (1) (1) (1.4) dl 
- 0".0 
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5 Cl Cl 

3a $ 




Ksin (Lap - 180°) C = SI sin a B secd B da B 

a 

(.00097) (1) (.02676) (1) (1.4) da. 

0 " .0 



Di d x 



= K COS ( Le P - ^ 



O to 

+ 180 ) sina B secd B 



3 Di d x 

31 



di 



3Pi di 
3 dm 



dd B 



- 0".0 



= K cos (Lap - n + 180°) sina m tand 0 secd n dd B 
= (.000097) (1) (1) (1) (1.4) dd. 



3 Did! 

3o\j 



3 to 

d^ m 



o ,f .o 

o Q 

= Ksin (Lap - to + 180 ) cos a B secd B da 

= (.000097) (1) da? 

= 0".0 



We wall now consider the standard error in lunar right ascension. We 
will use the quadratic form of equation 2.3, eliminating terms which are 
zero. We will assume that da? and 



3Bibi 

3P 



dp 



are dependent terms while all other non-zero terms are independent. 
Therefore, 



to 2 
(d a gapp ) 



+ ^ da ?dp 

dp 



and 



to 

d a sapp 



= 400 tan 2 d B cos 2 a? + 100 
+ 100 tan 2 d B cos 2 a? 

- 200 tan 2 d 0 cos 2 a? 

= \Zs00 tan 2 dmCos 2 a? + 100 
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Ai a[ 



d 

+ ( P 0 1 + a ) sin I cosa m 



da 
5 a 



a 

= sin I cos a B da 

a 



= . 02676 cosa m da 

Max. value = 0",5 (will be neglected) 



dAi ai 

5 1 



a 

dl = (P 0 t + a ) cos I cos a 0 dl 



= - .0015059 dl 

= 0".0 



dAi al 
Sd D 



dd m = 0".0 



SAi ax’ , Cl _ d _ Q O 

F) d 0! m = - ( P 0 1 +a) - sin I sina m da n 
5a V 



Q 

= (-.0015059) (-.02676) (1) (da m ) 



Bi bi 



a 

= p sina n 



— , Bl bl . d p= - sina ^ dp 
dP 



5 Bx b 



5 I 



-dl = 0".0 



8 d 1 .^ dd, = 0".0 



5a? 



a n 

Pcosa n da „ 
- .00029 da? 
0".0 



Ci cl = -Ksin(I^ P - £2 + 180°) cos I sina 

-K sin (Lq P - £2 + 180°) sinl cosdm 



£ 



sin d m 
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dl 



S Ci cl 



dl 



= ( + K sin ( Lq P - £2 + 180°) sinl sina a sin dm 
- K sin ( Lap - £2 + 1S0° ) cos I cos d = ) d I 
= -.00007 dl 
= 0 M .0 



3 Ci cl 
Sd a 



ddni 



( - K sin ( Lap - £2 + 180°) cos I since a cos dm 
+ K sin (Lap - £2 + 180°) sinl sind a ) dd a 
-.00007 dda 
0".0 



3 a? 



da 



£2 



- Ksin ( Lap - £2 t 180°) cos I cos a n sind m 
.00007 da? 

0".0 



Di dl 



api dl 

ai 



dl 



apT d* 

a dm 



ddn 



- K cos (Lap - £2 + 180°) cos a 0 sindm 

0".0 

- K cos ( La P - £2 + 180° ) cos a ? cos d 0 

.00007 ddm 

0".0 



aDxdi 1 

a a ? 




= Kcos (Lap - £2+ 180°) sina B sindm 
= .00007 da? 

= 0 M .0 
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We will now consider the standard error in lunar declination. We will 
use the quadratic form of the similar declination equation 2.3, eliminating 
terms which are zero. We will assume dd m and 





9 Bi bi d p 

dp 


are dependent terms. 




Therefore, 

( dd sapp ) 


... .s . .9 Bib! , .2 
- (ddm) + ( 9p dp) 

d Bibi' , „ u 
+ . dpdd B 

bp 

2 Cl s Cl 

- 400 sin ol a + 100 sin a m 

2 Cl 
-200 sin ol m 

2 Cl 

= 300 sin ol m 


and ddsapp 


Cl 

= 17". 3 sinQ! „ 



Thus, we can see that if stars with near zero declination are used, we 
can effectively limit the error in lunar right ascension. Let us examine 
the circumstances if we use stars with ± 10° declination and ± 10° right 
ascension. The maximum error will then be: 



J Cl 

d OL s a pp - 


fk maGJSMtsi III 1 WMH ' IMMlWIlWii 1 ' 1 

: V 300 (.03109) (1) + 100 = 10”5 


ddsapp ~ 


* 17". 3 (.17633) = 3". 1 



Thus, we can see that there can be an optimum time for which to obtain 
a position. If we consider limiting to values of 350° right to 010° 

and 170° right to 190°, we will immediately limit the error in declination 
to a maximum value of 3".l. See Appendix E for stars which fall into 
this area. 
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Next we will consider errors in time. We have defined 



LMDT + eq.N = LADT 
or 

RMDT + eq.N = RADT 

We can consider that RMDT is errorless so that the error then occurs 
only in eq.N which is dependent upon CT and T. So that we may say 

dT 3 = da 2 + dT 2 + dadT 

or dT = 10 s . 0 draconic time or 5".0 of arc. 

We will now consider the standard error of each method in deter- 
mination of astronomic position. 

For Method 1, we will first present the differential formula for 
latitude [16] : 

d4> = - sec A d Z - cos'&tanAdh 

- sin $ cosec Z sec d - tan d cot Z dd 
cos A 

We desire the azimuth to be near 0 or ISO degrees so that 

d 1 ? = - dZ - (sin^cosec Z sec d - tan d cot Z ) dd 

and if Z equals approximately 45°, then for a station at or near the equator 
d mil also equal about 45° so that 

d4> = - dZ - dd 

Therefore, for Method 1, these errors are independent, letting dZ = 7”.0, 
and considering 

(d$>) 2 = (7’’.0) 2 + (17". 3 sina?) 2 

and d$ = 7”.0 to 18". 7 

Thus, an error in h will have little effect on latitude determination. The 
error in latitude determination is basically related to the approximate right 
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ascension. The differential for azimuth [16] is 

dA = ( sin^cosec 3 Z - sin d cosec Z cot Z ) dh 
+ cos <I>sin h cosec 3 Z dd + cot Z sin Ad$ 

If we consider the zenith distance near 90° we see 

dA = (sin<i>- sind) dh + cos<i>sinhdd 

If the station is near the equator, will be close to 0 . Therefore, 

dA = - sin ddh + sin hdd 

Finally, we will select a pair of stars close to the meridian, preferably 
one with 90° north declination and one 90° south. Therefore, dA = 0 
for stations near the equator when we combine the pair. If we use only 
a sextant, we will have no means of accurately turning the angle between 
a star and landmark so that no permanent azimuth could be computed. 

The equation for longitude [16] is 

Cl 

X = GWp + h - (T + ATn) 

the apparent right ascension 
the hour angle as determined 
the local time on time keeper 
the chronometer correction 
Cl 

= da sapp + dh - dTn - dAT n 

= - (cosec A sec <£dZ + cot A sec <£d <£) 

- (tan^cosec Z sec d - tan d cot Z sec d) dd 
sin A 

The error in h may be eliminated by observing at an azimuth of 90° or 
270° and by observing east and west stars in pairs which are symmetrical 
to the meridian and at the same altitude. Considering tan$ is zero, the 
star pairs will eliminate the declination errors and zenith error terms. 



where a 

h 
T 

AT n 

Therefore dX 
and dh 
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So that 



. ft 

dX = da app - dT - d AT n 
d AT n has little effect on the moon for any ordinary chronometer, so that 

dX = daf?pp - d T 

Therefore ' d>? = (da Sapp f + (dT) 2 - dTda^ pp 

p 3 ft 

= 300 tan 2 d 0 cos o; a + 100 + 25. 0 - 
5. y^OO tan 2 d a cos' 2 & < ‘n + 100 
Limiting d n to 45°, we see 

dX = 8". 6 to 18”. 0 

For Method 2, we have already stated that the low precision of the 
AOT precludes further consideration. 

For Method 3, with impersonal operation and high equipment precision, 
we can consider that the only errors presented will be systematic in nature, 
and these specifically will be the systematic errors in star coordinates 
and time. Thus, we can consider that an astronomic position by the use 
of this method will yield errors near the equator (declination < 10°) 

d<±>= dd^pj, = 3".l (for optimum time) 
dX = da|? pp - dT 

dX 2 = (d & s ap P ) 2 +(dT) 2 - dT da^pp 
dX = 9". 1 

and dA = sin Aj cot Z d$+ cos § (tan - cosAj cotZ)dX 

near the equator, we say 

dA = (sin Aj d<±>- cos Aj dX) cot Z 

The largest Z we can hope for is 8°, and assigning an azimuth of the j th 
slit to be 45° , we find: 



37 



dA 2 = 25 d'S' 2 + 25 dX 2 (considering no dependence) 
Therefore, dA > 50". 0 

For Method 4, the errors are basically the same as for Method 1, 
if a panoramic camera is used, and basically the same as for Method 3, 
if the zenith camera is used (considering camera operating at optimum 
time). See Table 3 for comparison of results. 



Method 



d$ 



dX 



dA 



1 

2 



7". 0 - 18". 7 

> 20". 8 



8". 6 - 18". 0 

> 20". 8 



not directly 
obtainable 
not directly 
obtainable 



3 



3".l (optimum time) 9".l > 50". 0 



4* 



3". 1 - 18". 7 8". 6 - 18". 0 > 0".0 



♦Dependent upon horizon or zenith camera type, not including 
instrument or observation errors. 



Table 3. - Estimated Standard Errors (Station Near the Equator) 



2. 6 Conclusions 

Due to the possibility of limiting the error in star coordinates , it 
would seem that an automatic instrument as described in Method 3, would 
be the one which would be able to take best advantage of this point. 
Missions, of course cannot be scheduled so that these star coordinates 
could be utilized. It would be only coincidence if the current short 
duration landings occured at the proper time. Thus, the advantages of 
an instrument which could determine position automatically after the 
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astronauts depart is apparent. During future missions, if the landings are 
near the lunar equator, it would seem that the Electronic Zenith Camera/ 
Zenith photography would take best advantage of the stars with near zero 
declination. If landing were far removed from the equator, optimum 
times would need to be recomputed due to the increased error in. lunar 
right ascension, which would result from the higher star declinations. 

The format of the current missions, therefore, dictates that we use an 
automatic instrument which could take advantage of stars which offer 
minimum coordinate errors. The ease of set up and minimum training 
required for the astronauts, combined with the inherent accuracy of the 
system, show that the Electronic Zenith Camera or similar instrument is 
the best possible choice of equipment for the current missions. An 
accurate azimuth can only be determined by the use of an horizon camera 
at this time. 

It must again be stated that the Electronic Zenith Camera has not yet 
been manufactured, and could not be included in immediate missions. 
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3. DEFLECTION OF THE VERTICAL 



Deflections of the vertical relating the normals of the astronomic 
and selenographic systems can be obtained by: 

(1) Direct astronomical measurement of latitude and longitude 
with subsequent comparison with the selenodetic coordinates 
extended from a datum by triangulation or traverse. The 
results are affected by accumulation of error in the 
triangulation or ter verse. 

(2) Integrating gravity anomalies over large areas surrounding 
the point where the deflection is required. 

(3) Computing the deflections of the vertical due to topography, 
and then refining this system using some model of isostatic 
compensation. 

(4) Computing the deflections of the vertical from known spherical 
harmonics. 

The first method depends on an assigned deflection at a datum, thus 
is of no use in this paper. The second method requires gravity anomalies 
over large areas, and these of course are not available except by obtaining 
them through method four. But, method four is of little value in deflection 
computations as the harmonics reflect moon-wide rather than local 
distribution. The computation of deflections of the vertical requires more 
local data near the station. Thus, we will investigate method three, but 
will use method four as a check on the magnitude of the deflections. 

In addition, the effect of other bodies on the deflection of the vertical 
will be calculated. 

3.1 Computation of Topographic Deflections of the Vertical 

Following the method of Jordan-Eggert [11] adapted for the moon, we 
consider the influence of the visible masses in the neighborhood of a station. 



For this, we imagine the terrain around the station divided into compart- 
ments by means of concentric circles and radial lines, so that we obtain 
a great number of vertical columns between a given inference surface 
and the moon's physical surface. 

If the components of this deflection in the meridian and the prime 
vertical are expressed by A£ and Ap respectively, then we have 



3 0 H P ,T 

AC = 77 ; (sina 2 - sin«i) log e r 2 / n 

4 O' 7T ma 

- 3 q h p ,r 

( cos ot 2 - cos a i) loge r 2 / r i 

4 C> It B a. 



where © 



H 

a 

a x a z 
r i , r 2 



The surface layer density 

The mean density of the moon 

Height above reference surface 

The radius of the moon 

The compartment azimuths 

The outer and inner radius of compartments 



The reference surface which will be used is defined as a moon centered 
sphere of radius 1,737.988 kilometers. 

We will use the following values for the constants as given in [ 6 ]. 

0 = 2.670 gm/cm 3 (assumed)* 

© m = 3.343 gm/cm 3 

a = 1,737,988 meters 

p = 206,265" 



* This figure was obtained from averaging figures in various sources 
for a surface layer of 10KM. An error in this figure will produce 
a proportional error in the deflections. 
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and we obtain 



A£ = 0". 0226 H(sino: 2 - sin^i) log e r 2 / r x 

Ary = 0".0226 H( cos 0:3- cos a 1) log e r 2 / r 1 

Note: H in meters 

These two values are to be computed for each terrain section around 
the point 0, and then the two components of the deflection of the vertical 
for the point 0 are equal to the algebraic sum of all values of A£ or of 
A 77 as: 



We have so far assumed that the whole of the attracting stratum is 
in the horizon of the station, which for the moon would be grossly in 
error when only a short distance from the station. Therefore, it is 
necessary in order to secure the required degree of accuracy to take 
into account the curvature of the moon's surface. The effect of this 
curvature is to throw the compartment surface below the horizon of the 
station. Therefore, when significant, we will use the equation developed 
by A. R. Clarke , given in [11] . 



n 



n 



4 = s A£i 



T) = S Ar]! 



i-i 



A 4 = 0".0226 H (sin^s - sinOi! ) log e 



tan l/4 ® g 
tan 1/4 © 1 



+ cos 1/2 0 3 = cos 1/2 © 1 



A 77= - 0".0226 H (cos a s - cos a 1) log e 



tan 1/4 Q 2 
tan 1/4 © 1 



+ cos 1/2 © 3 - cos 1/2 © 1 



in which © 3 is r 2 reduced to an arc 
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radians 



© = r 2 ( in km ) 

2 1,737.988 

and ©i is ri reduced to an arc 



©i 



ri ( in km ) 
1,737.988 



radians 



Any values for either the constant r 2 / ri or for the constant difference 
(sinG! 3 - sin^i) may be adopted arbitrarily. Advantage is taken of this 
fact to adopt such values as would make the computation as simple and 
rapid as possible. The values r 2 / ri = 1.19362053 and (sinQ! 2 - sin<2i) = 
0.25 were adopted. With these values, the deflection equations will be: 



A£ = .001" h (in meters) 

At) = ,001"h (in meters) (the reference line is shifted 90°) 



As a result then, of this particular selection of arbitrary constants 
defining the limits of the compartments, the deflection produced at the 
station by the material lying above the reference surface in any compart- 
ment is expressed as one second of arc for each kilometer the mean 
elevation is above the reference surface. 

The arbitrary selection of the adopted values of (sino: 3 - sinQ! i ) 
and r 2 , ri was guided by three considerations. First, it was important 
to be as rapid as possible in determining the contribution of each section. 
The compartments must be small enough to bring the accuracy of the 
method within required limits and yet large enough to avoid excessive 
amounts of detail in the computation. The compartments should be 
compact areas, not long and narrow, in order to facilitate the estimation 
of the mean elevation within each compartment. 

The radii of the adopted circles separating the compartments are 
shown in Table 4. Each group of sixteen compartments (four in each 
quadrant) forms a ring. The rings, for convenience of designation, 
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were assigned serial numbers commencing with the inner ring, (see 
Figure 8 ) . 

It was determined that the curvature of the moon had to be taken 
into account at ring 46. The outer radii of all succeeding rings were 
computed iteratively on the IBM 360-75 computer, (see Appendix A). 

It was found by Hay ford [8] , that on the earth, the computed 
topographic deflections are much larger than the actual deflections of the 
vertical, in that an influence is in operation which produces an incomplete 
counterbalancing of the deflections produced by topography, leaving much 
smaller deflections in the same direction. The deflections of the vertical 
are due to irregularities in the distribution of the masses composing the 
earth. These irregularities may occur either as a result of irregularities 
in the surface of the earth (topography), or as a result of irregularities in 
the distribution of the densities beneath the surface. 

The distribution of the densities below the surface of the earth is 
invisible and largely unknown, but the theory of Isostasy postulates pre- 
cisely such a relation between sub-surface densities and surface elevations. 

If the earth were composed of homogeneous material, its figure of 
equilibrium, under the influence of gravity and its own rotation, would be 
an ellipsoid of revolution. 

The earth though is composed of heterogeneous material which varies 
considerably in density. If this heterogeneous material were so arranged 
that its density at any point simply depended upon the depth of that point 
below the surface, or, more accurately, if all the material lying at each 
equipotential surface (rotation considered) was of one density, a state of 
equilibrium would exist and there would be no tendency toward a rearrange- 
ment of masses. 

If the heterogeneous material composing the earth were not arranged 
in this manner at the outset, the stresses produced by gravity would tend 
to bring about such an arrangement, but because the material is not a 
perfect fluid, the rearrangement will be imperfect. In the partial 
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RING 


Outer Radius (KM) 


RING 


Outer Radius (KM) 


1 


.0080 


36 


3.9211 


2 


.0095 


37 


4.6803 


3 


.0114 


38 


5.5865 


4 


.0136 


39 


6.6682 


5 


.0162 


40 


7.9592 


6 


.0194 


41 


9.5003 


7 


.0231 


42 


11.3398 


8 


.0276 


43 


13.5354 


9 


.0330 


44 


16. 1562 


10 


.0393 


45 


19.2844 


11 


.0470 


46 


23.0184 


12 


.0561 


47 


27.4754 


13 


. 0669 


48 


32.7956 


14 


.0799 


49 


39. 1461 


15 


.0953 


50 


46.7267 


16 


.1138 


51 


55.7757 


17 


.1358 


52 


66.5781 


18 


.1621 


53 


79.4741 


19 


.1935 


54 


94.8707 


20 


.2310 


55 


113.2545 


21 


.2757 


56 


135.2082 


22 


.3291 


57 


161.4304 


23 


.3928 


58 


192.7602 


24 


.4688 


59 


230.2080 


25 


.5596 


60 


274.9951 


26 


.6680 


61 


328.6057 


27 


.7973 


62 


392.8567 


28 


.9517 


63 


469.9958 


29 


1. 1359 


64 


562.8462 


30 


1.3558 


65 


675.0283 


31 


1.6184 


66 


811.3203 


32 


1.9317 


67 


978.2927 


33 


2.3057 


68 


1185.5146 


34 


2.7522 


69 


1448.1099 


35 


3.2850 


70 


1793.1006 



Table 4. - Ring Outer Radius 
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FIGURE 8. 



Lunar Template 



rearrangement some stresses will still remain as different positions of 
the same horizontal stratum may have somewhat different densities, and 
the actual surface of the earth will depart from the ellipsoid of revolution 
in the sense that above each region of deficient density there will be a 
bulge or bump, and above each region of excessive density there will be a 
hollow. The bumps on the earth would be mountains and the hollows 
would be ocean floors. The excess material represented by that portion 
of the mountains above the reference surface will be compensated for by 
a defect of density in the underlying material and vice-versa for the 
oceans. This, of course, can be explained by the famous Pratt-Hayford 
theory of Isostasy. 

If we now consider the moon, we see that the consequences of 
hydrostatic theory as stated by Kopal [15] demand that given a sufficiently 
long time the material which on a short time-scale may behave as a 
solid is bound to get crushed under its own weight to settle to a form of 
minimum potential energy, which is a sphere. The departures of the 
lunar globe from spherical form permissible for a figure of equilibrium 
in the prevailing field of force are indeed small. To this extent, the 
moon may seem to be in hydrostatic equilibrium. However, this agree- 
ment becomes noticeably worse when we compare the momenta of the 
lunar globe with their theoretical values. Finally, by considering the 
mechanical ellipticity, we see that its value is so much at variance with 
the requirements of hydrostatic equilibrium that the discrepancy must be 
considei'ed as real. 

Therefore, in spite of the approximate prevalence of hydrostatic 
equilibrium in the lunar interior, as attested by the nearly spherical form 
of its globe, the moon's motion around its center of gravity reveals the 
presence of small but unmistakable departures from this equilibrium. 

On the other hand, O'Keefe [19] was able to show by comparison between 
calculated Bouguer anomalies and the essentially free air anomalies of 
Muller and Sjogren, that the topography of the moon could be isostatically 
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compensated and that to a rough first approximation the moon is in 
isostatic equilibrium. He considers an actual crustal layer on the 
moon, for which he states the Pratt-Hayford theory of Isostasy could 
be applied. 

Under the assumption of Pratt's hypothesis of Isostasy, we will 
introduce a surface of compensation at the depth T below the reference 
surface. If we denote the changes of A£ and A 77 due to the changes in 
topography by A £ 0 and At 7 0 , then the sums A£ + A £ 0 and Arj + Atj 0 
are the terms which, when summed, give the total deflection of the 
vertical. From Jordan- Eggert [11] : 



A£o 



- 3 A e T p" 
4 7 ! © B A 





(cos a 2 - cos a 1 .) log e 



r 2 + V T 2 + r 2 2 



n + V t ? + n 3 



A ©f = ©H 



Therefore: 




(sino: 2 Sinai) [log e — 




let F 




log, 
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Then for a definite value of T we can set up a table, from which for 
eveiy value ri (because r 2 = constant Xr! ) the factor F can be 
taken. Then we will have 

A£ + A £ 0 - FA£ 

At) + At) 0 = fAt) 

O'Keefe states [19] that for 20% difference in lunar crust and lunar 
mantle density, a depth of compensation of 25 kilometers is in order, 
while for 10% difference, a depth of compensation of 50 kilometers is 
in order. Further extending these assumptions, calculations for depths 
of compensation of 25, 50, 75 and 100 kilometers has been made using 
the program and subprogram in Appendix A. 

Table 5 gives the various F factors for each ring owing to a 
specific assumed depth of compensation. 

Calculations have also been made for 125, 250, 500, 1,000 and 1,500 
kilometers depth of compensation primarily to obtain a full set of values. 

There is direct evidence of excess mass in the maria over what 

would be called for by Isostasy. The explanation of this excess mass 

which was discovered by Muller and Sjogren of JPL is subject to various 

interpretations. O'Keefe [19] states that this excess mass may be due 

to the maria being filled with some sort of material which arrived after 

the formation of the highland topography. H. C. Urey [21] states that 

this mascon effect could be produced by a flat circular slab of chrondritic 

material. For instance, Urey calculated that in Mare Imbrium, this 

disc would be some 4 kilometers thick and be 670 kilometers in diameter 
2 1 

(mass = 5 X 10 grams). 

J. G. Stipe [20] states that these mas cons could be caused by meteorites 
which have bombarded the moon and have been of necessary size and at 
necessary depth to create these mass anomalies. 
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2 

3 

4 

5 

6 

7 

8 

9 

10 

11 

12 

13 

14 

15 

16 

17 

18 

19 

20 

21 

22 

23 

24 

25 

26 

27 

28 

29 

30 

31 

32 

33 

34 

35 



25 km 


50 km 


75 km 


100 km 


0.9997 


0.9998 


0.9999 


0.9999 


0.9996 


0.9998 


0.9999 


0.9999 


0.9995 


0.9998 


0.9998 


0.9999 


0.9994 


0.9997 


0.9998 


0.9999 


0.9993 


0.9996 


0.9998 


0.9998 


0.9992 


0.9996 


0.9997 


0.9998 


0.9990 


0.9995 


0.9997 


0.9998 


0.9988 


0.9994 


0.9996 


0.9997 


0.9986 


0.9993 


0.9995 


0.9996 


0.9983 


0.9991 


0.9994 


0.9996 


0.9979 


0.9990 


0.9993 


0.9995 


0.9975 


0.9988 


0.9992 


0.9994 


0.9971 


0.9985 


0.9990 


0.9993 


0.9965 


0.9983 


0.9988 


0.9991 


0.9958 


0.9979 


0.9986 


0.9990 


0.9950 


0.9975 


0.9983 


0.9988 


0.9941 


0.9970 


0.9980 


0.9985 


0.9929 


0.9965 


0.9976 


0.9982 


0.9915 


0.9958 


0.9972 


0.9979 


0.9899 


0.9949 


0.9966 


0.9975 


0.9879 


0.9940 


0.9960 


0.9970 


0.9856 


0.9928 


0.9952 


0.9964 


0.9828 


0.9914 


0.9943 


0.9957 


0.9795 


0.9897 


0.9932 


0.9949 


0.9755 


0.9878 


0.9918 


0.9939 


0.9708 


0.9854 


0.9903 


0.9927 


0.9651 


0.9826 


0.9884 


0.9913 


0.9584 


0.9792 


0.9861 


0.9896 


0.9504 


0.9752 


0.9834 


0.9876 


0.9408 


0.9704 


0.9802 


0.9852 


0.9294 


0.9646 


0.9764 


0.9823 


0.9158 


0.9578 


0.9718 


0.9789 


0.8996 


0.9496 


0.9664 


0.9748 


0.8804 


0.9399 


0.9599 


0.9699 



Table 5. - F Factors 
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lrn 

36 

37 

38 

39 

40 

41 

42 

43 

44 

45 

46 

47 

48 

49 

50 

51 

52 

53 

54 

55 

56 

57 

58 

59 

60 

61 

62 

63 

64 

65 

66 

67 

68 

69 



F Factors (Continued) 



25 km 


50 km 


75 km 


100 km 


0.8577 


0.9283 


0.9521 


0.9641 


0.8309 


0.9145 


0.9429 


0.9571 


0.7994 


0.8981 


0.9319 


0.9489 


0.7627 


0.8787 


0.9188 


0.9390 


0.7200 


0.8556 


0.9032 


0.9272 


0.6712 


0.8285 


0.8847 


0.9133 


0.6163 


0.7965 


0.8627 


0.8966 


0.5558 


0.7593 


0.8368 


0.8769 


0.4908 


0.7161 


0.8063 


0.8535 


0.4232 


0.6669 


0.7707 


0.8260 


0.3557 


0.6115 


0.7293 


0.7936 


0.2910 


0.5505 


0.6818 


0.7559 


0.2319 


0.4852 


0.6282 


0.7122 


0. 1803 


0.4175 


0.5687 


0.6624 


0.1370 


0.3501 


0.5045 


0.6065 


0.1022 


0.2858 


0.4373 


0.5451 


0.0751 


0.2272 


0.3695 


0.4795 


0.0545 


0.1763 


0.3040 


0.4118 


0.0392 


0.1338 


0.2435 


0.3445 


0.0280 


0.0996 


0.1902 


0.2806 


0.0199 


0.0731 


0.1452 


0.2225 


0.0141 


0.0530 


0.1086 


0. 1722 


0.0099 


0.0381 


0.0800 


0.1304 


0.0070 


0.0272 


0.0582 


0.0969 


0.0049 


0.0193 


0.0418 


0.0710 


0.0035 


0.0136 


0.0299 


0.0513 


0.0024 


0.0096 


0.0212 


0.0368 


0.0017 


0.0067 


0.0150 


0.0262 


0.0012 


0.0047 


0.0105 


0.0185 


0.0008 


0.0033 


0.0074 


0.0130 


0.0006 


0.0023 


0.0051 


0.0091 


0.0004 


0.0016 


0.0035 


0.0063 


0.0003 


0.0011 


0.0024 


0.0043 


0.0002 


0.0007 


0.0016 


0.0029 


0.0001 


0.0005 


0.0011 


0.0019 
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Wise and Yates [22] state that mascons are produced by mantle plugs 
upwelling into giant impact bases punched through the lunar crust, followed 
by volcanic filling of the remainder of the crater above the plug. 

No matter which theory is accepted, the direct result will be a further 
deflection of the vertical. The area where the mascon is located, will 
deflect the vertical, proportionately larger than the height of this area 
would indicate. The exact numerical increase in the deflection would be 
dependent upon the mass distribution as offered in the theories above. 

We could consider these mass increases as additional point masses and 
compute their effect on the deflections accordingly, or we could select a 
station which would minimize their effect. 

Let us consider a central station which is quite distant from the major 
mass anomalies so that their effect will be insignificant (very low F factor) 
in the overall result. For instance, considering any station near the 
equator from 20° east longitude westward to 40° west longitude will limit 
the contribution of these mass anomalies on the deflection of the veitical 
as the station will be greater than 600 kilometers from the outer edge of 
these anomalies. We can see for a station greater than 600 kilometers 
from a large mass anomaly that the largest F factor is only about . 015 
(down to 100 km compensation depth) and the contribution of the mass 
anomalies are then generally negligible on the overall result. For 
topographic deflections only, the area of these mass anomalies will 
generally fall into one compartment only, and will in turn affect the 
result very little. For this reason, I have selected the center of the 
crater Landsberg "A" as a starting point. Its coordinates are 328° 57’W 
and 0° 13'N. This point was also selected due to availability of charts. 

Due to the nearly concentric area in terms of height of the center of 
Landsberg "A” and the use of a 1:250,000 chart, I have neglected rings 
1 to 35, which are approximately the inner 3 kilometers of rings. Rings 
36 to 50 were computed using the AMS series topographic lunar map sheet 
1 edition 1 of 3 - 64. This map is a modified stereographic projection 



52 



with a point of origin 4,410.10 kilometers from the point of tangency of 
the map plan and a selenographic sphere with a radius of 1,737.988 
kilometers. The contour interval is 250 meters based on a vertical 
datum centered at Mosting "A" at an elevation of 7,000 meters. The 
horizontal datum was also based on the center of Mosting "A" which was 
considered to be latitude 3° 10’ 47" S and longitude 354° 50' 13". Rings 
51 to 63 were computed using an edition 1 - AMS sheet 2 topographic 
lunar map, scale 1:2,000,000 with the same datum, etc., as above. 

Rings 64 to 68 were computed using an edition 2 - AMS sheet 2 
topographic lunar map, scale 1:5,000,000 with the same datum, etc., 
as above. 

Rings beyond 68 were neglected due to their approach to the lunar 
limb; and their coinciding reduction in accuracy of height determination. 

Tables 6 and 7 give the heights of each compartment for computations 
in the prime vertical and in the meridian respectively. 

Table 8 shows the computed deflections for assumed depths of 
compensation. 
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Table 6. - Compartment Heights for Computations in the Meridian 
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Table 6. - (Continued) 
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Table 7. - Compartment Heights for Computations in the Prime Vertical 
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Table 7. - (Continued) 



Compensation Depth 
25 km 


+ 


2". 934 


V 

+ 10". 983 


50 


km 


+ 


3". 332 


+ 20”. 709 


75 


km 


+ 


2". 234 


+ 27". 757 


100 


km 


+ 


0". 883 


+ 33". 123 


125 


km 


- 


0". 547 


+ 37". 447 


250 


km 


- 


6". 118 


+ 49". 548 


500 


km 


- 


16". 792 


+ 1' 11". 984 


1,000 


km 


- 


30". 184 


+ 1* 38". 947 


1,500 


km 


- 


38". 529 


+ 1* 47". 427 


CD 




- 


1’04". 450 


+ 2’ 17". 480 



Table 8. - Deflections Considering Depths of Compensation 



We cannot consider that the value of © will remain constant down to 
the depths which we are considering. Thus, as depths increase, the 
average density for this surface layer will vary. We consider that the 
average density for these depths are as follows in Table 9 (interpolated 
from Kaula [12] ) . 
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KM 


Average Density 


FI 


25 


2.670 


1 


50 


3.000 


1.1236 


75 


3.140 


1.1760 


100 


3.210 


1.2022 


125 


3.250 


1.2171 


150 


3.263 


1.2221 


250 


3.336 


1.2494 


500 


3.343 


1.2521 


1,000 


3.300 


1.2360 


1,500 


3.343 


1.2521 



Table 9. - Average Surface Density for Selected Depth 



The FI factor in Table 9 is a deflection multiplier, and its effect 
on deflections is given in Table 10. 
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Compensation Depth 






V 




25 


km 


+ 2". 934 


+ 


10”. 


983 


50 


km 


+ 4”. 118 


+ 


23". 


269 


75 


km 


+ 2". 627 


+ 


32". 


642 


100 


km 


+ 1". 062 


+ 


39". 


820 


125 


1cm 


0". 666 


+ 


45". 


447 


250 


km 


7". 644 


+roi". 


905 


500 


km 


- 21". 025 


+ 1*30". 


131 


1,000 


km 


- 37". 307 


+ 2’02". 


298 


1,500 


km 


- 48”. 242 


+ 2*14". 


509 


00 




- 1'04".450 


+ 2’17". 


480 



Table 10. - Deflections Considering Changes in Density for 
Changes in Depth of Compensation 



The graph in Figure 10 shows the relationship of depth of compensation 
to deflections, considering no change in surface density, and considering a 
change in the surface density for changes in the depth of compensation. 
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Figure 9 . - Deflections Considering Changes in Density for Changes in Depth of Compensation 



The maximum topographic deflections on the moon will most likely 
occur at a highland -maria dividing line. For instance, a north-south 
deflection was examined at the southern wall of the crater Pitatus at 
latitude 31° 10’ S and longitude 346° 10' . With a few assumptions (due 
to lack of lai-ge scale charts on the area) it was determined that this 
deflection was approximately 5 minutes. It is quite possible that other 
larger* deflections exist, but probably will not greatly exceed this value. 

It is interesting to note that the deflections may be greatly minimized by 
choosing a station which is surrounded by nearly concentric rings of equal 
elevation. This will yield nearly zero deflections for the inner zones. 
This ideal station could be at the center of a concentric crater, at the 
top of an isolated rise on flat ground, or at a central point of a flat 
area. If, of course, a station were selected which did not satisfy this 
concentric condition, the inner rings must be considered. 

3.2 Computation of Spherical Harmonic Deflections of the Vertical 

In order to compare the values obtained by the method of isostatic- 
topographic deflections, we will use the formula as given by Dragg [3] to 
compute the deflections by the use of spherical harmonics. 

00 

jL N p ^ 

£ =— TT s s (AnmcosMX + Bnh sin M X) — — Pnm (sin (p) • 

ycoscp N= 2 [f=o LN - 1 

sincp - 

*7 = ^ 7 n S S TTT (Ann Sin M X - BnmCosMX) P nm (sin<p) 

/ COS KJ fjr2 JN *" 1 

Dragg [3] neglected the 20, 21, and 40 spherical harmonic terms in 
his computation for the earth, but they have been included in this paper. 



- M + 1) 
N - 1 



P N + i , * ( sincp) 



62 



I will use the spherical harmonics coefficients given by Blackshear [6], 
(see Appendix D). These are C n m and Snm coefficients which must be con- 
verted to A n m and B n m above. 

By definition we see that Cnm = Jnm and Sn M = Knm . 

For Heiskanen and Moritz [9], we see 

A*nm = GM r N Jnm 

B’ N m = - GM r N Knm 



and that from Dragg [3] 



* - SI La a A , 

A n m - , A m 

N + 2 

r 



Bnm — ~ ^ 7 ^ B’nm 



N t 2 



Thus we see that in terms of coefficients Cnm and Snm 



£ = £ £ (C n m cos MX + S N MSinMX) ](N + 1) P N M(sin<p) • 

COS<P n= 2M t 

sin <p - (N - M + 1) Pn + i« m (sin<p)] 

CO 

- 1 N 

7? = £ £ M(CNMSinMX - SnmCOsMA) P N M(sin<p) 

cos <p n =2 m=o 



This equation was evaluated for the Landsberg "A" station using the 
IBM 360-75 computer, the program of which is contained in Appendix B. 
The results of this evaluation are 

£ = - 12*. 184 
7? - 2'13". 707 

These values fall between the isostatic-topographic results of 250 km 
depth of compensation, and the results of topographic -compensation only, 
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when surface density changes are not considered. These values fall 
between the isostatic-topographic results of 250 lan and 1,500 km depths 
of compensation, when surface density changes are considered. 

3. 3 Effect of Other Bodies on the Deflection of the Vertical 

It is obvious that the attractive effect of the sun and the earth will 
affect the deflection of the vertical on the lunar surface. The following 
developments refer directly to the earth but apply equally as well to the 
sun. Planets also have an effect but for all practical purposes this effect 
can be considered insignificant. 

Following the method of Jordan - Eggert [11] adapted for the moon, we 
see that 



Deflection 



3 M e A 3 sin Z cos Z 
M m r 3 



where 



M e - The mass of the earth 



Mo - The mass of the moon 



A - The lunar radius 



Z - Zenith Distance of the earth 



r - Earth-moon distance 



and where 



= 81.301 
M, 



A' 



3 



(1737.63 f 
(384, 400) 3 



and 



s 

r 



9.27 X 10' 



Therefore: 

Deflection = 2.253 X lO^p" sin Z cos Z 



4". 646912274 sin Z cos Z 



2". 323 sin 2 Z 
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For the effect of the sun we will use the formula 



Deflection 



Where 



Ms 

M n 



and 



A 

“3 

r 



Deflection 



3 Ms A 3 sin Z cos Z 

— 3 

Mm r 



27,158,000 



( 1,737.63 f 

(149,415, 834. 4f 



1.5728 X 10' 



0". 026 sin Z cos Z 
0". 013 sin 2 Z 



which is negligible in comparison with other uncertainties introduced. 

We will now divide the deflection of the vertical caused by the earth’s 
attraction into a northward and westward directed component. If A is 
the azimuth of the earth counted from the north, then we have for the two 
components of the deflection of the vertical 

i e = 2". 32 sin 2 Z cos A 
r] e = 2". 32 sin 2 Z sin A 

These deflections are not constant for a station but vary in aceoi’danee 
with the change in earth-moon distance and the geometric librations in 
latitude and longitude. They also vary with the physical librations but 
these are insignificant. 



3.4 Conclusions 

Due to lack of knowledge of the lunar interior and its isostatic 
properties, we cannot, with any degree of confidence, state the magnitude 
of the deflections of the vertical at a given point. 

For the computation of the topographic deflections, we used lunar 
charts with a lunar centered spherical surface as the equipotential 
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reference surface, and from these the compartment relative heights were 
determined. If we consider the effect of an offset sphere on this computation, 
we can see immediately that the relative heights would change in accordance 
with the translation of the sphere's center. For instance, if the center 
of the reference sphere were offset toward the earth, the relative heights 
of the eastern compartments calculated in this paper would be consider- 
ably less, and the rj deflections, considerably smaller. Thus, we will 
need better height determination on the moon , before the topographic 
deflections can be considered to be sufficiently accurate. Due to the 
relative ease of computation of the spherical harmonic deflections, it may 
be best to use them at present instead of the topographic deflections, 
especially in areas which have local concentric heights. The computed 
deflections could be, of course, corrected for the effect of the earth. 

If in the future. Method (1) is employed, and relative deflections 
between stations were obtained, we would be in a much better position to 
predict the absolute deflections at these same stations using the topo- 
graphic method. 

The employment of Method (2) would, if extensive enough, yield the 
absolute deflections of each station. It is considered that in this method, 
lies the best possible results. But, alternately, the most time and effort 
would be required over the other methods. 

Further work can be accomplished using Method (3). For instance, 

Karl Jung in Report No. 41 of The Ohio State University, Institute of 
Geodesy, Photogrammetry and Cartography, has given the formulas to 
obtain approximate values of the influence of lunar craters on the 
deflection. These computations could become particularly useful for 
stations offset from the crater centers, as they would yield a rapid 
evaluation of the deflection attributed to the inner rings of the template. 
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APPENDIX A 



F FACTOR COMPUTATION PROGRAM AND ASSOCIATED 
MOON CURVATURE SUB - PROGRAM 

The main program computes the "F" factor for a given ring, given 
the depth of compensation. The associated sub-program computes the 
change in outer ring radius due to moon curvature. The main program list 
of variables is: 

IR Indexing variable for ring number 

J Indexing variable and integer compensation depth 

A1 Inner radius of ring 

A2 Outer radius of ring 

T Real number compensation depth 

FA Compensation depth squared plus outer radius squared 

FB Compensation depth squared plus inner radius squared 

FC Outer radius plus the square root of FA 

FD Inner radius plus the square root of FB 

FE FC divided by FD 

FF The natural log of FE 

FG The outer radius divided by the inner radius 
FH The natural log of FG 

F The F factor equal to 1 - FF/FH 

The sub-program Mooncu uses an iterative procedure employing the 
secant method to compute the outer ring radius taking moon curvature into 
account. 

The sub-program list of variables is: 

THETA 1 The inner radius divided by the lunar radius 

expressed in radians 

THETA 2 The outer radius uncorrected for lunar curvature 

expressed in radians 

THETB2 An approximate outer radius corrected for 
curvature expressed in radians 
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T1 

T2 

T2B 

T21B 

T21 

Til 

T3 

T3B 

SOC 



ASOC 



THETC2 



CSOC 

BSOC 



THETA 1/4 
THETA 2/4 
THETB2/4 
Tangent of T2B 
Tangent of T2 
Tangent of T1 
T21/T11 
T21B/T11 

The natural log of T3 minus the cosine of \ 
THETA 1 plus the cosine of ^THETA2 minus 
the natural log of the outer radius divided by 
the inner radius which is given by definition. 
This is considered the first approximation. 

The natural log of T3B minus the cosine of \ 
THETA 1 plus the cosine of §THETB2 minus the 
natural log of the outer radius divided by the 
inner radius which is given by definition. 

This is considered the second approximation. 

By the secant method where: 

x = Xt-T f(X,) - X, f ( Xl-T ) 

+1 f(Xi) -f(X 1 - 1 ) 



we see 



THETA2 (ASOC) - THETB2 (SOC) 

THETC2 = A»C 



The difference between X 1+1 and X! 
Outer radius in kilometers 
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A I \ 



GOC 1 DC 12 J = 25 f 12 b, 2 b 

0002 WRl Tt ( to, 70) J 

0003 7C FORMAT ( IrtCttoXt 'BASED CN~A OEWh OF " C JVREN S A T fO si uF',U, 

GOOA T = J 

0005 I ^ = 1 

0006 wRITt(6,*0) 

0007 80 FORMAT ( 9 X , 'RING* , 1 0 X , • F F AC TOR* , l C X , 'CUTER RADIUS*) 

0008 A 1 = . G 0 8 0 

0009 83 A 2 = I . lV3620'55*Al 

0010 3b hA=T**2+A2**2 

0011 FB= !**2 + Al**2 

0012 FC=A2+SQrU ( FA ) 

0013 F D = A l ♦ S Q R T ( F 5 ) 

00 1 A FE=FC/FD 

0015 " ' " F F = A LOG ( Fir ) 

0016 FG=A2/A1 

0017 FH=AlGG(FG) 

0018 F=1.-FF/FH 

0019 1 R = I R + 1 

0020 W <1 TF ( b ,dA ) I R , F , A2 

0021 d A r CRRA T { 9 X TTh , FOX , F 10. A , 3X f F 1C". Al 

00 22 A 1 - A2 

0023 IF ( 1R.LT.3C) OC TU 03 

002 A LALL ,VCC >CU( A1 , A2 ) 

0025 l F { IR.GT. 70) GC TO 12 

C026 GC TU 08 

0027 12 CONTINUE 

0C2 8 STOP 

0029 END 



f'C’CXCU 



ooci 




SUBROUTINE F'CCr-.CDt Ai ,-8SCC<l 


00C2 




f hr T A i = ~ 1 / 1 7 3 7 . & 3 


' 0003 


13 


r HF T A2 = 1 . 1 ) 36 20 5 3 * T n E T A 1 


OOOA 




THFTrf2 = 1.2-ThtTAl 


00C5 




I F t IrlL TAz.GT. 1 . ) GO TC "12 


0006 


1C 


T 2 = TH c T A2 / A • 


0007 




I 2S = ThFTo2/A. 


0003 




I UTnETAl/A. 


“0009 




T21R = IANTT2B) 


0010 




T21-TA\:(T2) 


0011 




T11=TA\(T1 ) 


0012 




I 3- I21/T11 


OC 1 3 




T3B = T21H/n 1 


00 1 A 




sac = a log t r 3) -cost 0 . eta 1 / 2 . 


ooT5 




A STTC = A’l L’ j ( r 3h ) - Cuj t THtT A l / 


0016 




THL TC 2 = ( I hr T A2- A SCC-T HE Tb2 


0017 




C SOC= T) E TC2-THET B2 


0018 




oSCC=l 737.63C*TnETC2 


0019 




I F ( AGS (CjOC) .LT. 1 .c-6) oC 


0020 




THE T A2 = Tj » t T 82 


“SO 21" 




T H t T ft 2 = T U E T~C 2 


0022 




GO IG 10 


0023 


1 2 


RMCR\ 


0024 




END 



) + CCS ( Thcl A2/2 . )-0. 1 7 
2 . ) FOS ( THE I 32/2. )-0. 
-SCO / ( ASOC-SUC ) 



TO 12 



699115 

1769911 



• K I LOME T • ) 
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APPENDIX B 



LEGENDRE POLYNOMIAL COMPUTATION PROGRAM AND 
SPHERICAL HARMONIC DEFLECTION PROGRAM 

The program computes the Legendre polynomials and the deflections 
of the vertical in both the meridian and the prime vertical using the 13th 
degree polynomial coefficients as given by Blackshear [6] for the position 
of the center of Landsberg "A". The Legendre polynomials were 
computed to the fourteenth degree as required and the subscript 15 
represents the Oth element [P(15,15) = P(0,0)]. The program list of 



variables is: 




P(15, 15) 


An array which contains the Legendre polynomials 


C(14, 14) 


An array which contains the Blackshear C 
coefficients 


S(14, 14) 


An array which contains the Blackshear S 
coefficients 


ALAT 


Latitude of computation point 


ALONG 


Longitude of computation point 


U 


Sine of latitude 


T 


Cosine of latitude 


I, J, KA, M 


Indexing variables 


MA,MB,N 


Indexing variables 


PSI 


Evaluation of a single term of the deflection 
in the meridian 


APSI 


Sum of PSI Terms 


PNU 


Evaluation of a single term of the deflection in 
the prime vertical 


APNU 


Sum of PNU terms 
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PPSI 

PPNU 



Multiplication of APSI by constants for final 
value of deflection in meridian 
Multiplication of APNU by constants for final 
value of deflection in prime vertical 
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000 i l) I Pc IS w:< <' ( 1 r i,l5l,C(K,H) ? S'(U,K) 

0002 _ !/t 2 0 > L 

000T K £ A U ( 5 , 2 U ) S 

0004 2C FUKF AT (4tlF.il) 

0 0 0 5 _ ALAT=.oG21£l6t>2 

0006 ALL) ,3 = -. 646416406 

0007 93 0 = 31 *( ALAT ) 

0008 r=CGS(_ALAT) 

“15009 ~ " J0'54 1=1,15“ 

0010 JO 63 J = 1 , 16 

0011 ?(I,J)=C.O 

0012 36 C0.\ f I »0E 

0013 64 l-OWT 1 \Ufc 

0014 P(15,l5)=I.O 

0015 >TTi 1 6 ) =0 

0016 P(2,16)=1.3*0**2— .5 

0017 P ( 2 , 1 ) =3. *0*1 

0018 P(2,2)=3.*T**2 

0019 UU 36 I=i,l4 

0020 P(I*15)=(( 2*1-1) *U*P ( I - 1 , l3)-( 1 - 1 ) * P ( I — 2 * 15))/I 

002 r p ( i , i T=C 2 *T-n*v*P( r-i , I — 1 ) 



0022 

0023 

0024 

0025 
C026 
002 7 
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56 



K A = 1 — 1 

P ( 1 , 1 ) ='P( 1-2, 1 ) ♦ ( 2*1-1 )*1 *P( 1-1 , 16) 
DU 67 J = 2 * K. A 

P( I , J ) = P ( I-2fJ) + C2*I-l)*T*P( I — 1 * J — 1 ) 

CONTINUE 

LOTT I NUE 



0028 


A PNU = 0 • 


0029 


A PS 1 =0 . 


0030 


JO 6 6 rj = 2 


0031 


UU b6 .V=l 


0032 * 


MA=K 


0033' 


Mb=M 


0034 


99 PS1=(J(M,' 




lni-K/.+ I ) *P 


0035 


PUJ=KA*P ( 


0036 


APS I = APS I 


0037 


APi4U=APNU 


0038' 


I F ( mYEC". 1- 


0039 


Fo CONTINUE 


0040 


H = 1 4 


0041 


ST A = 0 


0042 


UB=15 


0043 


GO TO 9 9 


0044 


bL’CONTi N!U'E~ 


0045 


PPS i = + AP;> 


0046 


PP \!U=- AP.>il 


0047 


V, 9 1 1 1 ( 6 , 3 i 


0048 


3 fi F0^i"AT(5>. 




1 T ICN Is Tl 


0049 " 


s Tup' 


(3050 


c NO 



V ) ^tS( ,«A*ALC,\lG) + 5 ( \‘ , .3 )*bl *( KA*AlUNG) ) * ( ( J+i > *P ( 4, MB ) *U- ( 

c n i , p> e ) ) 

T , r ’ b ) * ( C (N,V.)*S1 ••4(?-',A*ALCNG)-5(N,«)*COS(KA*ALaNG) ) 

♦PSI 
+ PfiU 

A ) GO Tu 85 



1*3600.0/(1*0.017453293) 
U*3oOO.C/ ( T*0. 01745 32 93) 



rut DcFLcC T 1G'\ I v) TriE fE^IClAN IS 
PjUHE 7 1 3 T I GAL IS * ,F*9.3///) 



* * F 9. 3/ / jX t 1 THE OEFLEC 
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APPENDIX C 



PHYSICAL LIBRA TION EQUATIONS 

As stated by Kolaczek [14] , physical libration is the physical 
deviation of the moon's motion from that described by Cassini's Laws. 

It can be treated as perturbations on the elements of the lunar orbit 
and rotation. Thus, the three components are called libration in 
longitude T , in node a , and in inclination p . The libration in 
longitude is a perturbation on the celestial longitude of the moon (can 
be considered as a correction to lunar time), libration in node on the 
celestial longitude of the node, and libration in inclination on the 
inclination I of the moon's equator. These perturbations are small, 
never exceeding 2' of arc. They can be expressed as follows: 

T = S at sin (ng + mg' + pF + qD) 

Icr = Ir +2bi sin(ng +mg' + pF +qD) 

p = Sc! cos (ng + mg' + pF + qD) 

Where g and g' are the mean anomaly of the moon and sun respec- 
tively (dependent on time). 

F = g + w 

D = g-g’+w-w' 

Where w and w' are the arguments of the perigrees of the moon and 
sun (dependent on time), n, m, p, and q and the coefficients a, b, and 
c are given in Table 11. 
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Symbols 


Series 


Argument Multipliers 
n m p q 


Coefficients in 
Arc Seconds 


T 


sin 


0 


0 


2 


-2 


a.1 — + 1.7 






0 


1 


0 


0 


+ 91.7 






1 


-1 


0 


-1 


- 1.2 






1 


0 


0 


-2 


+ 4.2 






1 


0 


0 


-1 


- 3.5 






1 


0 


0 


0 


- 16.9 






2 


-1 


0 


-2 


+ 1.0 






2 


0 


-2 


0 


+ 15.3 






2 


0 


0 


-2 


+ 10.0 


1(7 


sin 


0 


0 


2 


-2 


bi = - 3.2 






0 


0 


2 


0 


- 10.6 






1 


0 


-2 


0 


- 23.8 






1 


0 


0 


-2 


+ 2.5 






1 


0 


0 


0 


-100.7 


P 


cos 


0 


0 


2 


-2 


Ci = - 3.2 






0 


0 


2 


0 


- n.o 






1 


0 


-2 


0 


+ 23.9 






1 


0 


0 


-2 


- 1.9 






1 


0 


0 


0 


- 98.5 



Table 11. - Physical Libration Argument 
Multipliers and Coefficients 



77 



APPENDIX D 



BLACKSHEAR SPHERICAL HARMONIC COEFFICIENTS 

Table 12 below gives the C and S coefficients as determined by 
B lackshear in 1969. The coefficients are to the 13th order and degree. 
These coefficients are connected with satellite dynamics, with the 
potential V written in the form [9] : 

Km “ N a 

v = ~ {1 - S S (~) N [J nm R nm (©, X) +K nm S NH (0,X)] ) 

N= 1 M=0 



N 


M 


C 


S 


“ 0 


0 


1 .0000573551 £ + 00 


0. 


2 


0 


-2. 0IQ70 7 5 1 _66 E-_04_ 


o_._ 


3 


0 


-6 .30260 1_42 18 E- 0 6. _ 


0. 


4. 


0_ 


L.93_78_6_l 2.635 E- 0_5 


0. 


5 


0 


-7.4593773641 E-06 


o. i 


6 


0 


-1. 0784 697393 E- 06 


0. 


7 


0 


2. 40 33809155 E- 05 


0. 


8 


0 


— 2 . 6548898776 E- 05 


b. 


9 


0 


-1. 5429544143 E— 06 


0. 


10 


0 


5. 6344945 02 5 E- 05 


0. 


11 


0 


— 2.4602196863E-05 


0. 


12 


0 


3. 2989171 03 0E-05 


0. 


13 


0 


-5. 7723614632 c- 05 


o. 


2 


1 


-4. 4251 194184E-07 


-4. 5725447265E— 06 


3 


1 


2. 4370024 36 2 E- 05 _ 


1.. 2 . 30 0 6 37 6 7 07 E- 0_6 


4 


1 


-5. 732 38 04405 E- 06 


6.8017 765237 E-06 „ 


5 


1 


1 .36077442455-06 


- i . 2_I07 81 132 3 5-06^ _ 


6 


i 


-3.6076374435 E-0 7 


2 . 4 8 9 4 9 1 41 0 2 E - 06 


7 


1 


7. 82500 0287 6 E— 06 


8. 07 245 5 5222 E-06 


• 8 


1 


3.91 52106632 E-06 


- 1.929137C252E-06 


9 


1 


3. 11347120595-06 


- 5 . 33 30 3 865 4 7 E-06 


10" 


1 


9. 865666311 4E-07 


- 4 . 33 69 9077445-06 


11 


1 


-3.01 81 405386 E- 0 6_ 


. -.1.1.6 .2 29 1 3 9 6 9 E -05 


12 


1 


-7 . 49 i6 59 i 51 3 E-06 


1.741 194 3335E-C6 


13 


1 


— 1.2121962226 E— 0 6 


-5. 75S4594581E-06 



Table 12. - Blackshear Coefficients 
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6 6 _ - 2 . 0 2 6 7.4 6 3 32 6 E- 0 9 3. 89941632 12 E-0 9 

*7 6 - 1_. 0459807494 E-09 - 3 . 02 69 76 1 767E- 10~ 

8 6 3. 26 6789732 1 E-l 6 - 3 . 86978 8 275*9 E-jO^ 

9 _*_' 6 ljJJ 28 423 254 E-.’l 0 3_. 254*2 88 1 1 47E- 10 

"lb" 6 _ * 9_. 5 4 S2 2 1 5 37 0 E- 1 1 * " 6 . 06 5009*04*02 E^lY” 

jl _ 6 4 ._ 3 4 7 7 5J7 149 6 E- 1 1 _ - l . 71 88426 i98EMl_ 

. X 2_. . 6 . 4. 24. 8.2.6 J3_6_07 8 _E 3 1.1 3 . .7 £70 1 34 3. 9 8 E-_l 1_ 

13_ _ 6 4_. 367506896 8 £ 3 : 1_2 1. 533699 1203E-1 1_ 

7 _ _7 _ _1_. 20_2_6 8*1 677*9 ErLO_ ~3. 22 830*064 65JE- 10 

_jB_ ' __7 _ **_ 6 . 01 88 7 8*628 4 E-l 1 1 _ 3 .*59*32 703797 E— Il 

9 7 -_1. 802 35 50 33* 9_E 3 ’ii i . 0 5 6 1 0 8 1 7 8_3E - 1_0__ 

_1_0_ 7 _ -_i . 6 *1 2 ?_5 1 3 84 3 E- l 1_ _ - 2 .33 2 5 805068 E-J i_ 

_1 i 7 rT To 64.57 8083 E— 1 2_. . 3 1 . 07101 4 5.3±7E-1 1 _ 

12 _ _7 _ -1.54840332*73 E- 12* _ -2. 31 07720838E-1 2 

T~3~__~_7 _ 3 _ 1 _. 4 8 99 0 22.43 4 E - 1 2 2 . *4044*71 1 145E- 12 _ 

8 3 -1,. 7001 i 1699 4 E-l 1__ , . 3.00335333 7 9 E-l 1* 

9 8 2 ._ 0 5 84 3 5 G 2 7 5 E - 1 2 _ 7 ._88 3 474 3 2 0.7 EM 4_ 

10 *8 _1 . 5 8 9 .1.71 0 80 5 E- 1 2_ - 6 _. 5597619378E-12 _ 

J. \ 8 l _._01 I'Ll. IS 2.1 2 E- 1 _ 2 _ 8 .61 8.1 698 1*978- 13 

1 2 .3 1 _. 04 0071 91.86 £3 1.3 4.4 89ii03458E -i3 

*13 8 -_2 .JJ6 84 C 2_74 4 e'^ 14 5 . 57 1925*lj93 3Er 14_ 

9_ 9 9 ^ 61 7698**341 5 E- 13 - 1 .3*039 134543E r 12_ 

10 9 ..-1.9696887071 E.-.l 4. . -2 o 1.9 388.794 5 3E-13_ 

.11 9 - 1... 5 6.8 2.8.8 2 504 £-. 1 . 3__. 3. . . 5 . 2 ^ 2 01 36 2 6 S E-l 3. _ 

12 9. -7. 071205 7255 E - 14 - 3* . 55 3*4*89 C024E- 14 „ 



13 9 4. *68264698 74 E-l 6 _ z _l_.j92 2 j4 37 7 7J6E - 1_4 

10 __10 - 3_. 677 4006 812 E - 1*4 4. 65 57*2*0 1 1 26E- 14 

lT * _lTl_Ir 8TjT00 7jTq9 E- 1 5 J.oil29267p6E-i4_ 

12 *i 0_ 4 . 962 5 1*58_3 35 V- 1 5_ _ **- 1 *.*13 6 2 089*38 IE - 14 _ 

_1 3 1_0 _ 2_.Il|S0 r 2691*79E-T5_l " T. 472*261 0146E r l5. 

1 1 ll 3 .3 64*796 6 771 E- 1*6 ___-_l A 2 1J 23 53J38E- 15 _ 

12 11 7.4957 63 21 20 E— l 7 -2*. 4793959309~E~-16* 

13 *11 d -2.01 8 j~6 1 8 29 3 E- 1*6 * _1 . 7194186339E-16 

12 12* - 3.627 4.6 9 1 63 3E- 1 8 >_* - 1 . 5*5 262945 10E-18 

JL_3 _ _1 2 - 1. *87 62 16462 0 E-l 8 - 8_. 1 1 8902 1 267E- 18 

1 3*. 13* 9. 5 540605576 E-l 9*" "'* -1 .*70*5949*90 10E-19 .. 
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APPENDIX E 



OPTIMUM POSITIONS FOR MINIMUM 
STAR COORDINATE ERRORS 



Stars which yield the minimum error in transformation were found 
to possess the coordinates of: 

Lunar Declination = 0° 

o o 

Lunar Right Ascension = 0 or 180 

Thus this point will fall on the ecliptic at the position of the nodes. 
Figures 10 and 11 (from [16] ) give the positions of these stars in earth 
coordinates of right ascension and declination for the time January 0.5 
1971, January 0.5 1972, January 0.5 1973 and January 0.5 1974. The 
motion of these points is in direct relationship with the motion of the 
lunar nodes. 
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Figure 10. - Positions of Ascending Node 
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Figure 11. - Positions of Descending Node 
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